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We perform a model-independent analysis for recently reported data of the tt + tt invariant 
mass spectra in the T(45) — > T(15, 2S)tt + tt~ decays and point out that there does exist a 
broad peak below 0.6 GeV in the data of the T(45) — > T(lS)7r + 7r _ decay, which is analogous 



o 

to that in the T(3S f ) — > Y(1S)tt + tt~ decay. With the data of T(45) decays, we further test 



CD 



X 



our model developed for studying the puzzle in the T(3S) — > Y(IS)tt + tt decay. The result 



shows that with such a model, all the tt + tt invariant mass spectra of T(4S*) decays can 



be described. We also predict the cos#* distributions of T(45) — > T(1S, 2S)n + n decays, 
which can be used to justify our model prediction. 

PACS numbers: 13.25.Gv;12.39.Fe;12.39.Mk 

Keywords: T tt + tv~ transition, S wave tttv final state interaction, tetraquark 

T(45) is the first bottomonium state above the BB threshold. The branching fraction for 
T(45) — > BB is larger than 96% [lj. Recently, the 7r7r invariant mass spectrum of the T(4<S r ) — > 
Y(lS)ir + ir~ decay was measured by the Belle Collaboration [3, Q] and the BaBar Collaboration 
4], respectively, and the decay of T(4S) — > T(2S)ir + ir~ was reported by the BaBar Collaboration 
4|. Their results with rather large error bars are shown in Fig. [3j The branching ratios of decay 
modes reported by different collaborations are listed in Table U Although the branching ratio in 
Ref. [21 seems inconsistent with the newly reported one in Ref . _3|] , it is compatible with those in 



Refs 



An interesting phenomenon in the the Y(3S) — > T(lS')-7r + 7r decay is that the 7r + 7r invariant 



* Now at Institut fur Kernphysik (Theorie), Forschungszentrum Jiilich, D-52425 Jiilich, Germany 
^Electronic address: f.k.guo@fz-juclich.de1 



2 



TABLE I: The branching ratios of the decay modes T(45) -> T(15, 2S>+7r 





B(T(4S) T(lS , )7r+ 7 r")(10- 4 ) 


B(T(45) T(25)7r+7r-)(10- 4 ) 


CLEO [5] 


< 1.2 


< 3.9 


BaBar [4] 


0.90 ±0.15 


1.29 ±0.32 


Belle [2] 


1.1 ±0.2 ±0.4 




Belle [3] 


1.8 ±0.3 ±0.2 





mass spectrum has a double peak structure, namely a broad peak shows up below 0.6 GeV [(J. In 
past years, much attention has been attracted to such a phenomenon Now, 



the reported tt + tt invariant mass spectrum of the T(45) — ► T(2S)tt + tc decay also exhibits a 
e pea k structure, but in the T(4S I ) — > Y(1S)tt + it~ decay, the lower structure 



distinguishable doutx 

seems not evident [2j, 4, 14|. Based on this observation, Vogel suspected that the dipion transitions 
between two T(nS) states with An = 2 are special [151 ] . 

In order to explain the double peak structure and the angular distribution in the T(3S) — * 
T(1S)tt + tt^ decay, we proposed an additional sequential process in our model. In such a process, 
an intermediate bbqq state with J p = 1 + and 1 = 1 (called X) is suggested. By including this 
sequential process and considering the S wave tttt final state interaction (FSI), we well-described 
the T(nS) — > T(mS)Tr + ir~ (n = 2,3, m = 1,2 and n > m) data available in a systematical way 
Q]. Obviously, the newly reported data will provide us an opportunity to justify our model and 
check Vogel's assertion. 

In this paper, we perform a model-independent analysis for all the 7r + 7r~ invariant mass spec- 
trum data in the T(4S) — * T(1S, 2S)ir + ir~~ transitions. Although the Belle Collaboration claimed 

q n 

that the data in Ref. [2j should be superseded by those in Ref. [3(, we still use all these data 
to justify our model proposed in Ref. [7j, and consequently, predict the angular distribution for 
further confirmation in future. 

Before processing detailed calculations, we analyze the tt + tt~ invariant mass spectra of the T 
tttt transitions, qualitatively. 

1) The isoscalar character of bottomonia limits the isospin of the dipion system to be I = 0. 
Because both T(45) and T(IS') are the vector state, the tt + tt~ system favors the 5" wave than the 
higher ones. A mass difference of 1.12 GeV between T(4S) and T(15) gives the upper limit of the 
physical M n + n ~ region. In this region, the well-established isoscalar-scalar meson of /o(980) would 
couple to tttt strongly, which causes its width of about 40-100 MeV [l|. Therefore, there might be 
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a narrow peak or dip in the 7r + 7r~ invariant mass spectrum of the T(45) — * T(1S)tt + tt^ decay. 
However, if the bin in the experimental data is wide, as shown in all the reported data sets, the 
narrow state of /o(980) might not show up explicitly. 

2) In comparison with the measured invariant mass spectrum, namely the differential width, 
the squared decay amplitude (|.M| 2 ) reflects the decay dynamics more directly. Among low partial 
waves between two pions, only S- and D-waves can contribute to the decay amplitude. Because 
the contribution from the D-wave is suppressed by a factor of 2/45 due to the integration of 
the 2nd order Legendre polynomial, the squared decay amplitude can be approximated by the 
invariant mass spectrum divided by the integrated phase space. The ir + 7r~ invariant mass spectra 
of the T(35, 4S) — ► T(lS)7r + 7r~ decays without the phase space factor are shown in Fig. [U For 
comparison, we normalize the magnitudes around m w + w - = 0.41 GeV in all the data sets from 
Refs. 0,y,Q for the T(45) -> T(1S)tt + tt- decay and from Ref. for the T(35) -> T(15)vr+7r- 
decay to the values close to each other. From this figure, one sees that the tendency of the data 
in the T(45) — ► T(15)7r + 7r~ decay, especially in the region lower than 0.65 GeV, is analogous to 
that in the T(3S) — > T(1S)tt + tt~ decay. This indicates that in the tttt invariant mass spectrum of 
the T(45) — > T(15)7r + 7r _ decay, a broad peak does exist in the lower M n7T region, which implies 
that the above mentioned An = 2 rule might not be true. 



800 



600 



400 



200 



■I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

BaBar4S-1S 
+ Belle 4S-1S 
□ Belle (new) 4S-1S 
o CLEO 3S-1S exclusive 
• CLEO 3S-1S inclusive 



Normalized at -0.41 GeV 



1- 1 4: 



i .... i 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 

M . (GeV) 



FIG. 1: The ir + ir invariant mass spectra of the T(3S f , 4S) — > T(15')7r + 7r decays without the phase space 
factor. 



Now, we use our proposed model [7] to study T(45*) — > T(1S, 2S)7r + ir decays. The decay 
mechanism of T(nS) — * T(mS)ir + ir~ is shown in Fig. [21 where (a) and (c) represent the contact 



4 



diagram and the diagram in tree level with X, respectively, and (b) and (d) denote the correspond- 
ing diagrams with the S"-wave tttt FSI. 

i / / 

!/ + ,* + //+ 
^ u. / / i 

(a) (b) (c) (d) 

FIG. 2: The decay mechanism of T(nS) — > T(m5)7r + 7r _ . Solid external lines represent T(nS), dashed lines 
represent pions, and solid intermediate lines represent X. 



By incorporating appropriately the chiral expansion with the heavy quark expansion to the 
lowest order, the amplitude of the contact term in the heavy quarkonium ir + ir~ transition can be 



written as 



Vo = ~ T2 

J 7T 



(fllPi • P2 + 92PipI + 93ml)e* ■ e' + g^Pi^u + PiuP2^)e*^e" / , (1) 



where f n = 92.4 MeV is the pion decay constant, gi (i = 1,2,3,4) are coupling constants, p\ and 
P2 represent the 4-momenta of 7r + and 7r~, respectively, p® (i = 1,2) denote the energies of it^ in 
the lab frame, and e and e' describe the polarization vectors of the heavy quarkonia, respectively. 
It was shown in Refs. [ijlQ] that the g^-tevm in Eq. ([T]) can be ignored in the view point of QCD 
multipole expansion. Then, Eq. ([1]) can be reduced to 

y o = — iteiPi • P2 + g2p\p% + gsmDe* ■ e'. (2) 

J n 



Note that part of D-wave components still exist in the 52 term 



[71 



o„o _ 1 



P1P2 



1-/3 2 



pT 2 ~ Po(cos8;) - 2/3 2 pf P 2 (cos^ 



(3) 



where j3 is the velocity of the 7T7T system in the rest frame of the initial particle, p\^ = (p?*,Pi) 



and p\ = (P2*, P2) represent the four-momenta of 7r + and it in the cm. frame of the tttt system, 



in 



respectively. P (cos6>;) = 1 and P 2 (cos6£) = (cos 2 9* - l/3)/2 are the Legendre functions of the 
0-th order and the 2-nd order, respectively. 

In our model, we stressed that an additional sequential process T(nS) — > irX — ► T(mS)Trir 
should also contribute to the decay. Taking a simple S-wave coupling for the vertex T(nS) — > irX, 
the decay amplitude of the sequential process in Fig. [2] (c) can be written as 



V tree_ n . * I Zg^+^W^ , . ... 

VX — 9nm(-u e v — 2 7~- F r ~~2 2 ~- F~ ' V 4 J 

\p x + - m x + imxTx P x - ~ m x + im x i x ' 
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where g nm is the effective coupling constant among T(mS), T(nS), tt + and ir~ via an intermediate 
state X, and px± denote the momenta of the charged intermediate states of X^, respectively. 
In Ref. [7], we showed that the X state is necessary in getting a global fit to all the T(nS) — » 
T(mS)ir + ir~ (n = 2,3, m = 1,2, n > m) data. Especially in the T(3S) —* Y(lS')7r + 7r~ decay, 
the contribution of X is very important for both the dipion invariant mass spectrum and the 
angular distribution. The predicted mass and width of the X state are rax = 10.05 GeV and 
Tx = 0.688 GeV, respectively Q|. It should be mentioned that the predicted X state is supported 



171 ] and the QCD sum rules [I81 ]. 



by later theoretical calculations with the constituent quark model 
respectively. 

In Ref. [?|], it was also shown that the 5- wave -kit FSI is very important to reproduce the 
double peak structure in the 7r + 7r~ invariant mass spectrum of the T(35) — >Y(1S)tt + it~ decay. 
Such an importance has further been presented in the study of the heavy quarkonium chromo- 
polarizability that describes the interaction of heavy quarkonia with soft gluons [2^. The 
S*-wave 7T7T FSI can properly be treated by the so-called coupled-channel chiral unitary approach 



211 ]. In this approach, the tttt S- wave phase shifts can be well-described with barely one parameter, 
3-momentum cut-off q m ax = 1-03 GeV, and the low lying scalar mesons (a, /q(980), ao(980) and 



k) can dynamically be generated with reasonable masses and widths 
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22 1 . In the S- wave 



isoscalar sector, the tttt and KK channels are taken into account (for detailed information, refer to 



Ref. 



211] ). By using the phase convention \ir + ) = — 11,1) and the normalization (tt + tt \tt + tt 



(7r~7r + |7r~7r + ) = (7r°7r |7r 7r°) = 2 (considering the fact that 7r + ,7r~ and ir° are in the same isospin 
multiplet), the full amplitude of the it + -k~ + it~it + + 7r 7r° — > ir + ir~ process can be denoted by 
2t^j! c J r7r with t 1 ^ ^ being the full coupled-channel amplitude of the 1 = S-wscve tttt — > tttt process. 
In the T(4S) — > T(1S)tt + tt~ decay, the large phase space allows the KK channel to be on shell. 
Thus, the full amplitude of the K + K~ + K°K° — > it + it~ process can similarly be described by 
^kR tt-k/^^ with the phase convention \K° = — 1 1/2, 1/2) and the amplitude ^ which denotes 
the full coupled-channel amplitude of the 1 = 5-wave KK — > tttt process. To factorize the FSI 
from the contact term in Fig. [2](b), the on-shell approximation is adopted as usual. The off-shell 
effects can be absorbed into phenomenological coupling constants. In this approximation, both the 
tttt and KK loops shown in Fig. [2] (b) are allowed in the T(45*) — > T(15')7r + 7r~ decay, but only 
the 7T7T loop is permitted in the T(4S) — > T(2S)tt + tt~ decay. Then, the total transition amplitude 
for T(45) — > T(1S)tt + tt~ can be given by 

t = V + V 0S - G n ■ 24=^ + V 0S (m K ) ■ G 22 ■ 4^1U + V x™ + 9*^10% • 2t J =%^, (5) 
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and the transition amplitude for T(45) — > T(2S)ir + ir by 

t = V + Vos ■ G U ■ 2t T =% + Vt e + gne'^Gg ■ 2t{=%, (6) 

where Vos is the S'-wave projection of Vo, and Vqs(™>k) is the one for the process with two kaons. 
Ga and G^ are the two-meson loop integral {% = 1 for the tttt loop, and i = 2 for the KK loop) 
and the three-meson loop integral, respectively. 



^ (27r) 4 q 2 — mf + ie (p — p — q) 2 — mf + ie' ^ 



c ^_„. / A -r+PiPiK i ; (8) 

x 7 (2tt) 4 Fx ~~ m x + i e Q 2 ~ m w + * e (p ~ V ~ l) 2 ~ m2 v + 
where p' and p represent the momenta of T(45) and T(15), respectively. The loops are calculated 
with a cut-off momentum g ma x = 1-03 GeV which was used in explaining the data of the tttt S- wave 
scattering j^lj]. Thus, the FSI of tttt used here is consistent with the interaction of tttt employed in 
fitting the tttt scattering data. It should be mentioned that in calculating the transition amplitude 
for the diagram in Fig. [2Td), we do not consider the intermediate X s (bbss) state, the strange cousin 

n 

of X, due to its negligible contribution (refer to Ref. [7]). In this way, no additional parameter is 
introduced in the calculation. 

As has been pointed out in Ref. Q], to describe all the bottomonium dipion transition data, 
including the newly reported data, self-consistently, the 52 value is fixed to be 52/51 = —0.23. 
Then, the data of the tt + it~ invariant mass spectra in the T(4S') — ► T(15, 2S')7r + 7r - decays 0, [3], 
are fitted by three parameters, <?i, (73/51 and g nm , where m = 1 and m = 2 represent the final 
states of T(IS') and T(2S), respectively. Namely, in terms of the function minimization and 



error analysis package (MINUIT) 23J, these parameters can be determined by a least square fit. 
The resultant parameters with relevant errors are listed in Table HU In order to demonstrate the 
necessity of the suggested intermediate state X, the data fitting without X is also performed, 
and the relevant parameters are presented in Table [II] too. Moreover, the central values of the 
experimental branching widths, r(T(45) -» T(lS")vr+7r-) = 2.2 keV [2], 3.7 keV [3j, 1.8 keV [4] 
and r(T(45') — > T(25')7r + 7r _ ) = 2.7 keV are used in determining the physical values of coupling 
constants as well. It should be noted that due to g nm = g n x9mX, the values of g nm with different n 
and m are not fully independent. Using the values of 541 and 542 given in Table Ull and the value of 
531 given in Ref. [7J, one can deduce 532 = 7.36 GeV 2 which is different from the value of —0.00418 
GeV 2 given in Ref. Qj. Fortunately, because the contribution from the sequential process is not 
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important in the T(35) — ► T(25)7r + 7r _ decay, such a change does not affect the description of the 
data. 

The calculated tt + it~ invariant mass spectra of the T(4S) — > T(1S, 2S)it + ii~ decays are plotted 
in Fig. (3] where the solid and dashed curves denote the results with and without X, respectively. In 
the case with X, a clear bump exists in the low-energy region of the ir + Tr~ invariant mass spectra in 
both the T(45) — > Y(l«S)7r + 7r _ and the T(45) — > T(2S)tt + it~ decays, while in the case without X, 
no such a structure appears in the T(4S) — ► T(1S)tt + tt~ decay. It seems that the newly reported 
T(45) vr7r transition data support our model in which the sequential process with X is influential. 
Moreover, in the ir + ir~ invariant mass spectrum of the T(45) — * T(1S)tt + it~ decay, our model 
gives a narrow structure at the place just below 1 GeV. In addition, in the Y(4«S) — ► Y(2,S)7r + 7r _ 
decay, the calculated results for Fig. EJd) show a better description of the data if the sequential 
process with X is also considered. These results can be understood by analyzing the contributions 
of different terms. For simplicity, in the case of T(4S) — ► T(1S)tt + tt~ , we only focus our attention 
on the data set in Ref. [3|. The contributions from different terms are plotted in Fig. dj where 
the solid curve denotes the total result, and the dashed, dotted and dash-dotted curves represent 
the contributions from the terms without X, with X only, and the interference term, respectively. 
It is seen that the narrow peak close to 1 GeV comes from the S'-wave coupled-channel tttt FSI. 
Recalling that the isoscalar-scalar meson /o(980) can be generated dynamically in the coupled- 
channel chiral unitary approach [2l| , one can attribute the origin of the narrow peak to the effect 
of /o(980). However, due to the limited resolving power of the present data, this narrow peak 
might not be observed. Higher statistical data should be called. If this peak still cannot be found 
in the future higher statistical data, one should not be surprised, because it might be canceled by 
other stuff. An analogue is that in studying J/tp — > ujtt + tt~ decay in the coupled-channel chiral 



23], due to t 



24 



re contribution of the intermediate state &i(1235), the peak of 



251 ] . Similar to the situation in the T(35) — ► T(1S)tt + tt decay 



unitary approach 
/o(980) cannot be observed 
[7], the contribution from the sequential process plays a dominant role in the T(45) — > T(lS')7r + 7r _ 
decay. But, in the T(4S) — > Y(2S)ir + ir~ decay, the dominant contribution comes from the terms 
without X. Furthermore, in both T(4fi l ) decays, the contribution from the interference term is 
important. Anyway, in order to justify the heavy quarkonium tttt transition model further, the 
data with higher statistics should be requested. 

Using the central values of parameters in Table [Hi we calculate the cos 6* distributions for both 
T(45) dipion transitions (for Y(4«S) — > T(lS')7r + 7r _ , we use the parameter set determined from the 
data in Ref. [3J). The results are presented in Fig. [5] where the solid and dashed curves represent 
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TABLE II: Resultant parameters in fitting the ir + ir invariant mass spectra of the T(45) — > T(15, 2S)tt + tt 
decays. 





Decay 


Data source 


01 

.V 1 




Qi 1 Oi 
Ho/ i/1 


o nm (GeV 2 ) 






Belle [2] 


(6.00 ± 1.43) x 10" 


-3 


0.69 ±0.63 


5.28 ± 1.20 






Belle [2] 


(6.51 ± 1.18) x 10- 


-3 


0.44 ±0.48 


(fixed) 




-> T(lS , )vr+ 7 r- 


Belle [3] 


(5.56 ±4.48) x 10" 


-3 


0.94 ± 1.74 


7.60 ± 1.19 


T(45) 








Belle [3] 


(6.76 ± 1.46) x 10- 


-3 


1.04 ± 0.52 


(fixed) 






BaBar [4] 


(7.05 ±0.86) x 10- 


-3 


-0.94 ±0.53 


3.70 ± 1.02 






BaBar [4] 


(3.84 ±0.73) x 10- 


-3 


1.65 ±0.70 


(fixed) 


T(4S*) 




BaBar [4] 
BaBar [4] 


0.15 ±0.01 
0.14 ±0.01 




-3.37±0.17 
-3.67 ±0.14 


11.9 ±3.3 
(fixed) 




0.8 1.0 1.2 

M(tiV) (GeV) 




0.6 0.8 1.0 

M(juV) (GeV) 



1.2 



250 




0.8 1.0 1.2 

M(tcV) (GeV) 
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FIG. 3: The ir + ir invariant mass spectra for the decays T(4S f ) — > T(1S, 2S)tt + tt . The solid and dashed 
curves denote the results with and without X, respectively. 
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cose* cose* 



FIG. 5: Predictions of the cos6>* distributions, (a) T(45) -> T(lS>+7r-, (b) T(45) -> T{2S)ir+n-. The 
solid and dashed curves represent the results with and without the intermediate X state, respectively. 

the results with and without X, respectively. From the figure, one sees that the predicted cos 8* 
distribution of the T(4S*) — * T(lS')7r + 7r _ decay is much natter in the case without X than in the 
case with X. Clearly, the angular distribution of the T(4S) — ► Y(lS)ir + ir~ decay with X exhibits 
a similar behavior shown in the T(35) — > T(lS')7r + 7r _ decay, because the sequential process plays 
dominant role in both decays. Therefore, we can expect that the future angular distribution data 
would be a criterion for further judgment. 

Moreover, we would mention that in the T(4S) — > Y(15)7r + 7r~ decay, the numerical result 
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shows a negligible effect from the off-shell tensor resonance /2(1270). 

In summary, we analyze recently reported data of the ir + ir~ invariant mass spectra of the 
T(4S I ) — > T(1S,2S)tt + tt~ decays in a model-independent way, and find that the behavior of the 
7T7T invariant mass spectrum in the T(45*) — * Y(lS)ir + ir~ decay is an analogue of that in the 
T(3S) — > T(lS')7r + 7r _ decay, namely there does exist a broad structure in the tttt mass region 
below 0.6 GeV. This disagrees with the presumed An = 2 rule in Ref. [lj|. Then, we try to 
justify our model developed in studying the puzzle in the T(35) — > Y(IS)-k + tt~ decay by using the 
newly reported data of T(4S) dipion decays. We find that with the additional contribution from 
the sequential process proposed in our model, one can describe the experimental data, namely 
the intermediate state X plays an important role in the T(4S) — ► Y(1S)tt + tt~ decay. We also 
predict the cos#* distributions of Y(4S) — ► T(15, 2S')7r + 7r~ decays. The result shows that the 
cos#* distribution of the T(45) — * Y(1S)tt + it~ decay in the without X case is much flatter than 
that in the with X case. Therefore, this observable can be employed to justify our model in the 
future. 
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